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We consider the magnetic structure on the Fe(001) surface and theoretically study the scanning 
tunneling spectroscopy using a spin-polarized tip (SP-STM). We show that minority-spin surface 
states induce a strong bias dependence of the tunneling differential conductance which largely de- 
pends on the orientation of the magnetization in the SP-STM tip relative to the easy magnetization 
axis in the Fe(001) surface. We propose to use this effect in order to determine the spin character of 
the Fe(001) surface states. This technique can be applied also to other magnetic surfaces in which 
surface states are observed. 
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I. INTRODUCTION 

The properties of magnetic surfaces and interfaces have 
attracted recently a lot of attention because of the ad- 
vent of spintronics, a technology aiming to harness elec- 
tron's spin in data storage and processing, typically by 
utilizing heterostructures composed of magnetic and non- 
magnetic materials^. Recent theoretical studies revealed 
that the electronic properties of the Fe(001) surface 
can play a much more important role for ferromagnet- 
scmiconductor base d spintronic devices than it was pre- 
viously thoughtPsH. Using first-principles electron trans- 
port methods it was shown that Fe 3d minority-spin sur- 
face (interface) states are responsible for at least two im- 
portant effects for spin electronics. First, they can pro- 
duce a sizable Tunneling Anisotropic Magnetoresistance 
(TAMR) in magnetic tunnel junctions with a single Fe 
electrode^. The effect is driven by a Rashba shift of the 
resonant surface band when the magnetization changes 
direction. This can introduce a new class of spintronic 
devices, namely, Tunneling Magnetoresistance junctions 
with a single ferromagnetic electrode^ that can function 
at room temperatures. Second, in Fe/GaAs(001) mag- 
netic tunnel junctions minority-spin interface states pro- 
duce a strong dependence of the tunneling current spin- 
polarization on applied electrical bias. A dramatic sign 
reversal within a voltage range of just a few tenths of 
an eV was predicted. This explains the observed sign re- 
versal of spin-polarization in recent experiments of elec- 
trical spin injection in Fe/GaAs(001^P and related re- 
versal of tunneling magnetoresistcance through vertical 
Fe/GaAs/Fe trilayers. The TAMR effect was also ob- 
served recently in a Fe/GaAs/Au tunnel junction 7 . 

Many of the theoretical results mentioned above, are 
based on theoretical predictions that the Fe(001) sur- 
face band is of minority spin. However, to this date a 
direct experimental determination of the spin character 
of the Fe(001) surface band has yet to be done. Scan- 
ning tunneling microscopy (STM) is a well established 



technique for imaging surface structure d 8 * 9 * 10 * 11 ^ and for 
spectroscopic measureme nts of lo cal structures and in- 
homogeneities on surfaced- 12 * 13 * 14 !. Using this technique 
Stroscio et al provided the first experimental evidence 
that a surface band near the Fermi energy exists in the 
Fe(001) surface^, however the technique used was inca- 
pable of distinguishing the spin character of this band. 
Recent developments of STM tools with spin-polarized 
tip, SP- STM, a llow for controllable measurements of 
magnetid 16 * 17 * 18 ! and spin-dynamicaP^ features and lo- 
cal magnetic structured^. Spin-polarized tunneling for 
STM set-up has been theoretically studied with respect 
to noised, and spin-detection and spin- reversal of local 
spins located on a substrate surface^i 2 ^. 



In this paper, we consider a SP-STM set-up on top 
of a planar Fe surface in order to study the influence of 
the minority-spin resonant state at the surface on the 
spin-polarized tunneling current. We found that a com- 
plicated angle and energy dependence of the tunneling 
differential conductance emerges as a result of the en- 
ergy and momentum dependence of minority-spin band 
structure in Fe(001) surface. Our results identify a spe- 
cific route on how to determine the spin character of the 
Fe(001) surface band with the help of SP-STM. In prin- 
ciple, these results are applicable to a broad set of mate- 
rials where the minority /majority spin-structure exhibits 
nontrivial energy dependence. 



The rest of this article is organized as follows. In sec- 
tion II we discuss our numerical approach and present 
our results. Particular emphasis is given to the relation 
of our approach to Tersoff and Hamman formula for the 
tunneling current in STM^S. Then, based on the ob- 
tained results we discuss our proposal on how to detect 
the spin character of surface states with the help of SP- 
STM. Section III is the conclusion. 



FIG. 1: (a) Spin-resolved DOS for the bulk Fe (b) Spin- 
resolved DOS for the Fe(OOl) surface (c) Spin-resolved ky- 
integrated transmission for the Fe/vaccum/Cu as a function 
of energy. The Fermi level is at zero energy. 

II. MAIN 

We consider a Fe/vaccum/Cu tunnel structure with 
a nonmagnetic bcc Cu electrode. This electrode has a 
spin-independent free-electron-like band structure and a 
featureless surface transmission functional which makes 
it insensitive to the transverse wavevector. This elec- 
trode simulates an ideal STM tip. The structure in- 
vestigated consists of a semi-infinite Fe region, several 
layers of vacuum (empty atomic spheres), and a semi- 
infinite Cu region. Another advantage of such setup is 
that it avoids overall the possible 'handshake' of surface 
resonances at two opposite metallic surfaces/interfaces^. 
When two resonant states on opposite contacts are lo- 
cated in the same [E, kii ) space then naturally a resonant 
transmission equal to 1 will occur across the structure. 
Obviously a situation like this may occur in a symmet- 
ric structure. However, in real structures such 'hand- 
shakes' are unlikely because the symmetry of the struc- 
ture is broken by applied bias or disorder. Therefore, 
a Fe(001)/vacuum/Fe(001) setup may yield unphysical 
high transmission for surface states. 

The calculation approach is based on the Green's 
function representation of the Tight-Binding Linear 
Muffin-Tin Orbital (TB-LMTO) method in the Atomic 
Sphere Approximation (ASAJ2S. We use third order 
parametrization for the Green's function^. The elec- 
tronic structure problem is solved within the scalar rel- 
ativistic Density Functional Theory (DFT) where the 
exchange and correlation potential is treated in the 
Local Spin Density Approximation (LSDA). The con- 
ductance is calculated with the principal-layer Green's 
function techniqu e 1 28 * 29 * 30 ! within the Landauer-Biittikcr 
approach^. The semi-infinite Fe and Cu electrodes are 
separated by approximately 1 nm of vacuum represented 
by 6 monolayers of empty spheres. The structure is ori- 
ented in the [001] direction. Self-consistent charge distri- 
bution is achieved before any transport calculations are 



FIG. 2: Minority-spin ky -resolved DOS of the Fe(001) surface 
for different energies around the Fermi level. The abscissa is 
along [100] and the ordinate is along [010]. The maximum 
value is represented by red color, the minimum by blue. 

attempted. The spin-dependent ky -integrated transmis- 
sion 

T ff (E) = 1/2tt / t cr (E,k||)d 2 k|| (1) 

J 2DBZ 

is calculated in the window from Ep to Ep+eV. t CT (E, kii ) 
is the transmission coefficient and a =]\l ( T =majority- 
spin, j=minority-spin). The spin quantization axis lies 
along (001) direction. A uniform 250x250 mesh was used 
for the integration in the two-dimensional Brillouin zone 
(2DBZ). With this transmission we can associate a cur- 
rent density 

l-E F +eV 

J a (V) = e/h / T a {E)dE (2) 

Je f 

, this is an excellent approximation appropriate for com- 
parison with experiment when the considered voltages 
are small. Then the spin resolved differential conduc- 
tance is dJ a /dV oc T(E F + eV). 

Figs, [ija) and[TJb) show the calculated Fe bulk spin- 
resolved density of states (DOS) and Local DOS of the 
Fe monolayer at the Fe(001) surface, respectively. The 
energies are given with respect to the Fermi level Ep. 
In the Fe bulk (Fig. [TJa)) the majority-spin dominates 
over the minority-spin throughout the entire energy in- 
terval shown here. However, in the surface monolayer 
(Fig. [ijb)) the spin polarization of the DOS is totally 
reversed; it is the minority- spin that dominates over the 
majority-spin throughout the entire energy interval. In 
Refs. [5] and |3]it was shown that this reversal is caused 
by Fe 3d surface states of minority-spin. As shown in 
Fig. [ljc) , sign reversal of the spin polarization of the 
surface DOS doesn't necessarily lead to sign reversal of 
the spin polarization of tunneling transmission, at least 
throughout the same energy interval. 

To explain this in Fig. [2] we present the minority-spin 
ky-resolved DOS of the Fe(001) surface for different en- 
ergies around the Fermi level. The bright red features 



FIG. 3: The energy dependence of the differential conduc- 
tance of a FM tip as we rotate the direction of the magneti- 
zation in the tip by an angle 9 relatively to the magnetization 
in Fe. We show four different cases of DOS spin polarization 
P of the tip. The bright red regions observed for 8 = n are 
created by the minority spin transmission peaks in Fig. 1(c). 



on these plots are created by the surface band. The 
surface band has C^ v symmetry which is the symme- 
try of the Fe(OOl) surface. These bands are dominated 
by the minority-spin surface states arising from d x 2_ y 2 
and d xy orbitals on surface Fe sites that couple with the 
bulk Fe A2' minority band. Unlike the majority-spin 
bulk band this band never crosses the 2DBZ at the T 
point. The closer it gets to the T point is at the energy 
of Ep — 0.025ey where we see a bright four-petal struc- 
ture centered at the T point without touching it. Con- 
sidering that kj| is conserved during tunneling across an 
ideal surface, for an electron entering the vacuum region 
with a real k|| the decay rate of the electron wave func- 

whcre k is 



tion is proportional to exp 

the decay rate for normal incidence which is determined 
by the potential height of the tunneling barrier. The 
tunneling transmission for a given energy is also propor- 
tional to the total number of states at this energy n a (E). 
In Fig. [ljc) one can compare the spin-resolved tunnel- 
ing transmission with the spin-resolved surface DOS. We 
see the influence of both factors just mentioned above 
in the polarization of the transmission coefficient. The 
minority-spin transmission dominates over the majority- 
spin for a large part of the energy interval due to its much 
higher surface DOS. For the energy E = Ep — 0.025^, 
the minority-spin transmission has a maximum and in 
the energy interval where the minority-spin DOS is flat, 
the minority-spin transmission is less when the surface 
states are further away from the T point. For energies 
where the minority-spin surface states are far away from 
the center of the 2DBZ, the minority-spin transmission is 
less than the majority-spin even though for the same en- 
ergy the minority-spin surface DOS is much larger than 



FIG. 4: This is a three dimensional version of Fig3a) (P=l) 



the majority-spin. For bigger distances between the Cu 
counterelectrode and the Fe surface the ratio of minority- 
spin to majority-spin transmission should become even 
less for all energies in the shown interval. However, our 
calculations show that for a distance twice as big, the 
change of spin polarization bias dependence is not very 
big while the current drops by about four orders of mag- 
nitude. This shows that within the range of distances 
available to an SP-STM tip, the inversion of spin polar- 
ization should be detectable (below we describe how). 

We connect our results from the planar calculation to 
the realistic experimental STM set-up by the following 
observation. Following Tersoff and HamanrP^, the tun- 
neling current at K is given by 

I(V) cxj^f (HP)) [1 - / - eV)] x 

P k 



: |(i )P |T|i,k)| 2 5(^(p)-^(k)) 



(3) 



where / is the Fermi function and |(£, p|T|j,k)| is the 
tunneling matrix element between initial and final states. 
Here, we designate the momentum index p (k) to the Cu 
tip (Fe substrate) electrons and let i (j) denote the initial 
(final) state. We take the convention that electrons flow 
from the tip to the Fe surface for positive voltage V. In 
our planar case, the role of the tip is played by the Cu 
surface which, with regards to the tunneling, behaves as 
an STM tip. Namely, the electronic structure of Cu is 
known to be relatively featureless around Ep, therefore 
the tunneling matrix element in Eq. ^ only weakly de- 
pends on p. The substrate surface, on the other hand, 
has a strong in-plane dependence, while the out-of-plane 
component can be written as k = ky +q z . By integrating 
out the p dependence we can write (|3| as 



I(V) (x E k [1 
JdEjdk {l [1- 



- e y)]/(^(k))|(i|T|j,k)| 2 
eV)]/(^(k))|(z|T|j,k)| 2 



oc 
(4) 
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where in the coefficient of proportionality we have incor- 
porated the density of states of the tip (Cu) at Ep, and 
the kn integral is taken in the 2DBZ of Fc surface. Equa- 
tion ffl is equivalent to Eq. p| , while for a fixed energy 
Eq. HF is equivalent to Eq. (ITT. 

Based on these results, we would like to discuss the 
possibility of detecting the spin-character of Fe surface 
states with a ferromagnetic (FM) tip. Let's assume that 
the spin quantization axis in the FM tip is rotated by 
an angle 9 relatively to the spin quantization axis in the 
Fe surface. The spin components of the electron wave 



where, are the spin components of the transmis- 

sion coefficient (differential conductance) presented in 
Fig. [ljc) . When 9 = the spin quantization axis in the 
tip is parallel to majority-spins in Fe and when 9 = 180 
it is parallel to minority-spin. 

In Fig. [3] we show the energy dependence of the differ- 
ential conductance of a FM tip as we rotate the direction 
of the magnetization in the tip by an angle 9 relatively to 
the magnetization in the Fe. Four different cases of DOS 
spin polarization P of the tip are shown. In the general 
case, < P < 1, a complicated energy-angular depen- 
dence is observed. However, the trends can be easily 
understood from the two limiting cases of P = 1 (half- 
metallic tip) and P = (not shown in Fig. [3] but can be 
extrapolated from the case of P = 0.25). For P = 1 the 
energy dependence of the transmission for 9 = is iden- 
tical to the majority-spin transmission in Fig.Qc), while 
for 9 = 180 the energy dependence is identical to that of 
minority-spin in Fig. [IJc) (to facilitate the comparison in 
Fig. [4] we provide a three dimensional version of Fig. [3^). 
As we decrease the degree of spin-polarization in the tip 
the angular dependence of the differential conductance 
becomes less pronounced. As it should be, at the limit of 
P = the differential conductance becomes independent 
of the angle and equal to the sum + TK This case 
corresponds to the conventional (non-magnetic tip) STM 
measurement by Stroscio et al. 
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function in the tip can be written as 

I T, 6) = cos (9/2) I T) - isin (0/2) | |) 

I 1,0) =-i sin (9/2)\ T> + cos (9/2) \ [) (5) 

where, \u,9) (a =|,|) are the spinors for an arbitrary 
direction of the spin quantization axis and \a) are the 
spinors for the same spin quantization axis as in Fe. The 
total transmission coefficient for an arbitrary spin polar- 
ization P of the tip can be written as 



(6) 

I 

III. CONCLUSION 

In conclusion we have shown that unlike a conventional 
STM which only can measure a sharp peak in the energy 
dependence of the differential conductance when a local- 
ized surface band is present, an SP-STM measurement 
should be able to measure an angular dependence as well. 
In the general case of < P < 1 a complicated energy- 
angular dependence emerges but the trends can be easily 
understood from the limit of a half- metallic, P = 1, tip. 
The energy dependence of the differential conductance 
in this case will be monotonic for either parallel or anti- 
parallel direction of the magnetization of the STP tip 
relative to the direction of the easy axis in the Fe(001) 
surface. This can be used to extract the spin character 
of the surface band. 
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